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Abstract 



In this work we extend Varadhan's construction of the Edwards 
polymer model to the case of fractional Brownian motions in M. d , for 
any dimension d > 2, with arbitrary Hurst parameters H < 1/d. 
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1 Introduction 

In recent years the fractional Brownian motion has become an object of 
intense study due to its special properties, such as short/long range de- 
pendence and self-similarity, leading to proper and natural applications in 
different fields. In particular, the specific properties of fractional Brownian 
motion paths have been used e.g. in the modelling of polymers. For the self- 
intersection properties of sample paths see e.g. |GRV03] . |HN05] , |HN07j . 
[HNS08], |Ros87] . and for the intersection properties with other independent 
fractional Brownian motion see e.g. |NUL07j . [OSS10J and references therein. 
Comments on the relevance of fractional Brownian motion for polymer mod- 
elling, in particular with H = 1/3 for polymers in a compact or collapsed 
phase, can e.g. be found in [B C95] . 

The fractional Brownian motion on M. d , d > 1, with Hurst parameter 
H G (0, 1) is a d- dimensional centered Gaussian process B H = {B^ 1 : t > 0} 
with covariance function 

E(B^B^) = S f (t™ + s™ - \t - sn , i,j = l,...,d,s,t> 0. 

An informal but suggestive definition of self-intersection local time of a frac- 
tional Brownian motion B H is given in terms of an integral over a Dirac 
(^-function 

dt I ds5(B H (t)-B H (s)), 
o Jo 

intended to measure the amount of time the process spends intersecting itself 

in a time interval [0, T\. A rigorous definition may be given by approximating 

the (^-function by the heat kernel 

1 M 2 

(2Tce) d ' 2 
which leads to the approximated self-intersection local time 

L £ := f dt f dsp £ (B H (t)-B H (s)). (1) 

Jo Jo 

The main problem is then the removal of the approximation, that is, e \ 0. 

In the classic Brownian motion case (H = 1/2), L £ converges in L 2 only 

for d — 1. To ensure the existence of a limiting process for higher dimensions 



one must center the approximated self-intersection 

L £ , C :=L £ -E(L £ ). (2) 

For the case of the planar Brownian motion this is sufficient to ensure the 
L 2 -convergence of (J2J) as e tends to zero [Var69j, but for d > 3 a further 
multiplicative renormalization r(e) is required to yield a limiting process, 
now as a limit in law of 

r(e) [L e - E(L.)) . (3) 

Through a different approximation, this has been shown in |CY87j . |Yor85j . 

Extending Varadhan's results to the planar fractional Brownian motion, 
Rosen in |Ros87j shows that, for 1/2 < H < 3/4, the centered approximated 
self-intersection local time converges in L 2 as e tends to zero. 

This result, as well as all the above quoted ones for the classic Brow- 
nian motion, have been extended by Hu and Nualart in [HN05] to any d- 
dimensional fractional Brownian motion with H < 3/4. More precisely, Hu 
and Nualart have shown that for H < \jd the approximated self-intersection 
local time ([T]) always converges in L 2 . For 1/d < H < 3/(2d), a L 2 - 
convergence result still holds, but now for the centered approximated self- 
intersection local time fl2]). In this case, 

( TC H4 e- d / 2+1 l {2H ^ + o(e), if 1/d < H < 3/ (2d) 

E(4)= ' , (4) 

{ 2H { ly/i ln (±/£) + o(e), UH = l/d 

where Cn,d is a positive constant which depends of H and d. In particular, 
for 1/d < H < min{3/(2<f), 2/(d + 1)}, an explicit integral representation for 
the mean square limiting process L c as an Ito integral is even obtained in 
[HNS08J. For 3/ (2d) < H < 3/4, a multiplicative renormalization factor r(e) 
is required in |HN05] to prove the convergence in distribution of the random 
variable (J3J) to a normal law as e tends to zero. 

To model polymers by Brownian paths Edwards [Edw65] proposed to 
suppress self-intersections by a factor 

exp (-gL) , 

with g > 0. For planar Brownian motion Varadhan |Var69j showed that the 
expectation value E(L S ) has a logarithmic divergence but after its subtraction 
the centered L £ c converges in L 2 , with a suitable rate of convergence. From 



this, Varadhan could conclude the integrability of exp(— gL c ), thus giving a 
proper meaning to the Edwards model. For more details see also |Sim74j . In 
the three-dimensional case this is clearly much more difficult [Bol93j, [Wcs80j. 
In this note we extend Varadhan's construction to arbitrary spatial di- 
mension d > 2 and Hurst parameters H < 1/d. For this, the convergence 
results proved in |HN05j will be essential. Because of this, in the following 
section we collect from [HN05] the necessary information on fractional Brow- 
nian motion and its self-intersection local time, and in Section 3 we state and 
prove the existence theorem (Theorem [2]). 

2 Preliminaries 

As shown in [HN05] . given a <i-dimensional fractional Brownian motion B H 
with Hurst parameter H G (0, 1), for each e > the approximated self- 
intersection local time (JT|) is a square integrable random variable with 



(2n) d J T ((A + eXp + e)-/^/ 2 ' 



where 



T := {0, t, s', t') : < s < t < T, < s' < t' < T} 
and for each r = (s, t, s', t') G T, 

A(r) := (t - s) 2H , p(r) := (f - s'f H , (5) 

and 

/j,(t) := — [\s — t | + \s — t\ — \t — 1 1 — \s — s | J . (bj 

Furthermore, for each e, 7 > is 

E(L £ L 7 ) - E(L £ )E(L 7 ) =-- (7) 

1 



dr — : — - - — — - =: E, 



(2n) d J T \((\ + £ )( p + 1 )-^y/2 ^ X + £ )(p + 1 )y/2j ■ ^- 

Note that the integral in (JS} is also well-defined for all e, 7 > (however it 
might be infinite). Hence, using this integral representation, we can extend 
E £1 to general £,7 > 0. This is contrast with (JTj) which in general is not 
well-defined for 5 = and/or 7 = 0. 



From (jSJ) one can easily derive that a necessary and sufficient condition 
for convergence of L £C = L £ — E(L e ) to a limiting process L c in L 2 as e \ 
is that £"00 < 00. As shown in |HN05l Lemma 11], the integral Eqq is finite 
if and only if dH < 3/2. 

3 Theorems and Proofs 

Theorem 1 Assume that (d+l)H < 3/2, d > 2. T/ien t/iere exists a positive 
constant K such that 

E ((L £ , c - L c ) 2 ) < tfe 1 / 2 
/or a// e > 0. 
Proof. Using (jHJ), a simple calculation and taking the limit 7 \ yields 

E ((L e , c - L c f) = (E ££ - £ e0 ) + (£00 - £ £ o) 
with 

d 



Err — E, 



ee — ^eO 



2(27r) d 

1 



/ dr (A + e) I dx 



<0. 



!«A - :)(,>+. c)f 2+ ' {{\ + e){p + x)-^) d/2+] 
Hence 

E ((L £ , c - L c ) 2 ) < Eqo - £ e0 

= 2^?Jr dTP L ^ [(5 + xp)*/™ ~ ((A + x) P r +1 j ' ^ 

where 5 := \p — p 2 . Thus it is sufficient to establish a suitable upper bound 
for ([9]). Technically, this will follow closely the proof of Lemma 11 in [HN05J, 
based on the decomposition of the region T into three subregions 

Tn{s<s'} = T 1 uT 2 uT 3 , 



where 






{(t,s,t',s') 
{{t,s,t',s') 

{(t,s,t',s') 



< s < s' < t < if < T}, 
< s < s' < if < t < T}, 
< s < t < s < if < T}. 



Each substitution of T in OH]) by a subregion %, i — 1, 2, 3, yields a different 
case and for each particular case we will then establish a suitable upper 
bound. 

As in [HN05] , we will denote by k a generic positive constant which may 
be different from one expression to another one. We set D :— d + 1. 

Subregion 7i: We do the change of variables a := s' — s, b := t — s', and 
c = if — t for (t,s,t',s') G 7i- Thus, on 7i, for the functions A, p, and p 
defined in (Q and (Q we have 

A(£, s, £', s') =: Ai(o, 6, c) = (a + 6) 2H , p(£, s, *', s') =: pi(o, 6, c) = (6 + c) 2H 
//(*, s, £', s') =: m(a, b, c) = l - [{a + b + c) 2H + 6 2H - c 2 " - a 2H ] . 

On the region 7i one can bound ([9]) by the first term only, and to estimate 
the latter we shall use Lemma [3] below, yielding 

01 / £ rfx 1 <Ac 1/2 o 1/2 8' D/2 

Pl L (6 1 + x Pl )(^^- A£ Pl dl ■ 

From JHN05I eq. (59)], 

Sx > k(a + b) H (b + c) H a H c H > k(abc) 4H/3 , 
we deduce 

[ dadbdc5~ D/2 <k [ da db dc (abc)~ 2DH /3 <oo, 

because DH < 3/2. In conclusion the part of ([H]) stemming from integration 
over 7i is of order e 1 ' 2 . 

On the subregions 71, i — 2, 3, we have to consider the difference 



\ '■— Pi I dx I f! | _ ur)4 _ lW9 ,,, N jn-nirt I ' e ' : * " • 



(^ + ^) (zm)/2 ((A 4 + x) Pl ) (D+1)/2 



•m 



Subregion %: In this case we do the change of variables a := s' — s, 
b := i! — s' , and c = t — t' for (£, s, t', s') G 7i- That is, on % we will have 

X(t, s, t', s') =: A 2 (a, 6, c) = 5 2H , p(t, s, t', s') =: P2(o, 6, c) = (a + b + c) 
Mi, s, t', s') =: // 2 (a, b,c)= l - [{b + c) 2H + (a + b) 2H - c 2H - a 2H ] . 

In this case we decompose the corresponding integral (Q over the regions 
{b > rja}, {b > rjc}, and {b < rja, b < rjc}, for some fixed but arbitrary rj > 0. 
We have by ( TT6l) . see Appendix, 



/ dadbdcHf, < Cs 1 ' 2 \ dadbdcp 2 (A2P2) 

Jb>na Jb>na 



_ -j, 2 / dadbdc 



lb ^ a{a + b + c) DH b DH- 

If DH < 1, the integral is finite. If 1 < DH < 3/2, then by Young inequality 
dadbdcZf, < ke 1 / 2 f f dad ^ DH f dbb- DH 

b>T)a JO Jo [ a + c ) J-qa 

■T r-T 



< ks 1 ' 2 [ daa~ iDH / z+1 [ dcc- 2DH / 3 < 00. 
Jo Jo 

In the case DH = 1 we have 

/ dadbdcE £ 2 <ke 1/2 dcc~ 2/3 daa~ 1/3 ln(T/(rja)) < 00. 

Jb>rja Jo Jo 

The case b > rjc can be treated analogously 

To handle the case b < rja and b < rjc we first observe that 

for sufficiently small rj > 0. Hence, together with (Tl5l) . see Appendix, we 



obtain 



/ 



dadbdcE £ 2 < Ce 1/2 j da db dc p l 2 ' 2 p\ (A 2 p 2 ) _(Z?+2)/2 

< ke 1 ' 2 [ dadbdc (a AH ~ 2 + c m ' 2 ) (a + b + c) -™-™ b 2-2H-DH 

J b<ria,b<ric 



' b<ria,b<rjc 

< ke 1 ' 2 ! dadbdcb- DH (a + b + C y 2H - DH 

J b<r)a,b<rjc 
x / (2-£>/3)tf & Dff/3 + c (2-Z)/3)// & DH/3\ 

< ke 1 ' 2 [ dadbdcb~ DH (a + b + c)- 2H - DH a^ 2 - D ^ H b DH / 3 

< kE 1 ' 2 [ da db dc b -*DH/Z c -2DH,3 a -2DH/3 < ^ 

because DH < 3/2. 

Subregion 73: We do the change of variables a := t — s, b := s' — t, and 
c — t' — s' for (t, s, t', s') G % . Thus, on %, we have 

A(£, s, t' ', s') =: A 3 (a, 6, c) = a 2H , p(t, s, £', s') =: ^3(0, 6, c) = c 2H 

/i(t, s, t', s') =: ^3(0, b, c) = i [(a + 6 + c) 2 " + 6 2H - (6 + c) 2 " - (a + 6) 2 "] 

In this case we decompose the corresponding integral ([9]) over the regions 
{a > rjib, c > r/ 2 &}, {a < 771&, c < r/ 2 &}, {a > rjib, c < r] 2 b}, and {a < rjib, c > 
i] 2 b} for some fixed but arbitrary 7/1,772 > 0. By symmetry it suffices to 
consider the first three regions. Using ffl6l) . see Appendix, we obtain 



./ 



dadbdcE e 3 < Ce 1 ^ 2 / dadbdcp 3 (A 3 p 3 ) ' 2 

1 J a>mb,c>n9.b 



a>rjib,c>ri2b J a>rj\b,c>ri2b 

< ke" 2 I" db [ T ^f T 4^< ke" 2 f T J^< 00. 



n DH I r DH — I U2DH-2 

Jrub a Jr n b C JO u 



For the region {a < rjib, c < i] 2 b}, we observe that since H < 3/ (2D) < 1/2, 
we can conclude from (Tl5l) . see Appendix, together with [HN051 eq. (55)], 
i.e., /X3 < kb 2H ~ 2 ac, that 



% < Ce" 2 p l 3 /2 ^ 3 (X 3 p 3 ) 



-{D+2)/2 

< ke" 2 b iB ~ A a 2 ' ?u - DH c 2 - 2H - DH < ke" 2 a- 2DH ^c- 2DH ^b- 2DH /\ 



which is integrable. Finally, we consider the case {a > r]ib,c < 7726} • For 
772 > small enough we have 

"-5( (a+6r (( 1 + ^) 2H - 1 )- 62 "(( 1 + ^- 

< k ((a + b) 2H - 1 + b 2H - 1 ) c = kb 2H - 1 f(l + ^} 2H ~ l + l\c< kb 2H ~ l c } 

where in the last estimate we used IB. — 1 < (due to H < 1/2). Then using 
(|15p we obtain 

I dadbdc~l<Ce 1/2 I dadbdcpl /2 fj,l(\ 3 p 3 y {D+2)/2 

J a>r]ib,c<r]2b J a>rjib,c<r}2b 

< ke 1 ' 2 / da db b m ~ 2 a- 2H - DH / dc C 2 ~ 2H - DH 

Ja>rjib JO 

< ke 1 ' 2 f daa~ 2H - DH f ' ^ db b~ DH + 2H + 1 < ke 1 ' 2 f daa- 2DH+2 , 
Jo Jo Jo 

which is finite because DH < 3/2. ■ 

Theorem 2 (i) Assume that dH = 1, d > 2. Then there exists a positive 
constant M such that for all < g < M 

eM-gLc) (10) 

is an integrable function. 

(ii) Assume that dH < 1, d > 2. Then there exists 

L := limL e in L 2 

e\0 

and for all non-negative constants g 

exp(-gL) 

is an integrable function. 

Proof, (i) The case d = 2 and H = 1/2 was treated in [Var69j . In all 
other cases we are in the situation of Theorem [TJ In these cases we have a 

9 



logarithmic divergence of E(L e ) as e \ 0, see 01]). Combining this moderate 
divergence with the rate of convergence provided in Theorem [TJ the proof for 
integrability of ffTOj) for small enough non-negative g follows very close along 
the lines of |Var69[ proof of Step 3]. More precisely, by @, for < e < 1 
there exists a positive constant k such that 

T 

L e c > -E(L e ) > - k ; — -ttt I ln(e) I . 

For any constant N > k + ^hJ^kW^ I m ( e )l one nas 

P(^ c <-iV) = P(L C - L £ . c < -N - L e , c ) 

< Pf \L £C -L C \ > N- k ; — -— lln(e) 

V £ ' c c| - 2H(2n) d / 21 v ' 

An application of Chebyshev's inequality then yields 

P(L C < -N) < E " L "- L '' 2 » - 2 < K ^ , 

N - k -mk^ WE) \) { N ~ k ~ 5w*l ln ( £ >l) 

In particular, for 

e = exp (-H(2n) d/2 (N - k)/T) 

one obtains 

AK 
P(L C < -N) < ^3^ ex P {-H(^) d/2 (N - k)/(2T)) . 

Hence, there exists a positive constant M such that ( TTU1) is integrable for all 
< g < M. 

(ii) In the cases dH < 1 we know from |HN05j Theorem 1 (i)] that the 
following limit exists 

< L :=limL, in L 2 . 



e 



e\0 

Thus, exp(— gL) is integrable for all non-negative g. 



Appendix 

The following lemma is an immediate consequence of the Cauchy-Schwartz 
inequality. 

10 



Lemma 3 Let < a, (3 < oo and 1/2 < m < oo ; then there exists a positive 
constant A such that 

dx {a + (3x)~ m < Ae 1/2 a~ m+1/2 p- 1/2 . 
o 

For i = 2, 3 we set 

1 1 

The following lemma is a generalization of estimates (56) and (57) obtained 
in [HN051 Lemma 10]. 

Lemma 4 For i = 2, 3 i/iere exists a positive constant B such that 

Ux)<B^((\ t + x)p i )- {D+1)/2 - 1 , (11) 

ax)<B{{\ l + x)p i y (D+l)l \ (12) 

for all x > 0. 

Proof. Estimate (fiTj) implies estimate ( TT2l) . Indeed, according to |Hu01t 
Lemma 3 (2)], for some suitable constant < k < 1, 



Aipj - /ij = <5; > ZcAiPi. 

Since Aj,pj are positive, this implies that 

A < (1 - fe)Ai/0i < (1 - /c)(Ai + x) Pi , (13) 

for all x > 0. Thus, assuming ( 1TTT) . ( TT2|) follows from ( TT3l) . 
Therefore, the proof amounts to prove ( 1TT1) . Given 

^-(O-p^J -l|((A. + ,)«)-<^ (14) 

observe that due to ffTBl 

(Aj + a;)p 



0< 7—^ < i . -/•• < I. 



11 



Hence let us consider the function 

[0, 1 - k] 3 y H- f(y) := (1 - y)-( D+1 )/ 2 - 1 G [0, oo). 

Since /(0) = and /' is continuous on (0, 1 — k) with a continuous continu- 
ation to [0, 1 — k], there exists a positive constant B such that 

f(y)< max \f(z)\y< By for all ye [0,1 -A;]. 

«e[o,i-fc] 

Applying this inequality to f fT4l) yields the required estimate (ITT]) . ■ 

Lemma 5 For z = 2, 3 i/iere exists a positive constant C such that 

El < Ce^pTA (V,)- (D+2)/2 , (15) 

E!<Ce^ P ¥ 2 (\ iPi )- D/ \ (16) 

for all e > 0. 

Proof. Recall that 

5- = pi / dx£i(x), i = 2,3. 



Hence (ITBT) and (ITS]) follow from (ITT]) and (IT2]) . respectively, together with 
Lemma [U ■ 
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